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Nuclear reactions using rare isotopes has opened up several new frontiers in nuclear 
sciences [0.|3j]. In particular, heavy rare isotopes currently available at several radioactive 
beam facilities in the world and the more energetic ones to be available at the planned 
Rare Isotope Accelerator (RIA) provide a unique opportunity to explore novel properties 
of dense neutron-rich matter that was not in reach in terrestrial laboratories before. This 
exploration will reveal crucial information about the equation of state (EOS) of neutron-rich 
matter @|||j|]7]||] . Within the parabolic approximation (see e.g., |J), the EOS of neutron- 
rich matter of isospin asymmetry 5 = (p n — p p ) / (p n + p P ) can be written as 

e(p, S) = e(p, 0) + E sym {p)5 2 + 0{5 4 ) (1) 


where e(p, 0) is the energy per nucleon in isospin symmetric nuclear matter. The density- 
dependent nuclear symmetry energy E sym (p), especially at high densities, is very uncertain 
and has been a subject of extensive research with various microscopic and/or phenomeno¬ 
logical models over the last few decades. Besides being important for studying the structure 
of radioactive nuclei near neutron and proton driplines, the density dependence of nuclear 
symmetry energy has many profound consequences for several key issues in astrophysics 
0.0000. It is therefore an important theoretical challenge to make reliable predic¬ 
tions on which experimental observables can be used as sensitive probes of the E sym (p), thus 
constraining the isospin dependence of the nuclear EOS. Within a hadronic transport model 
it was recently shown that the tt~/ n + ratio is a particularly useful probe for the high density 
behavior of the E sym (p ) [|15|]. For central collisions of 132 Sn + 124 Sn a significant sensitivity 
of the 7r~/ir + to the variation of E sym (p) was found. How does this sensitivity depend on 
the isospin asymmetry of the reaction system? In this Brief Report, we investigate this 
question. An answer to it is not only interesting in its own right but also useful for planning 


future experiments 0. 

It is well known that the 7r~/n + ratio in heavy-ion collisions depends strongly on the 
isospin asymmetry of the reaction system, see, e.g., (0,0,00. It is also qualitatively easy 
to understand why this dependence can be used to extract crucial information about the EOS 
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of neutron-rich matter. On one hand, within the A resonance model for pion production 
from first-chance independent nucleon-nucleon collisions [^(J. the primordial tt~/ n + ratio is 
(5 N 2 + NZ)/(bZ 2 + NZ ) « ( N/Z ) 2 . It is thus a direct measure of the isospin asymmetry 
(. N/Z)dense of the dense matter in the participant region of heavy-ion collisions. It was shown 
that the ( N/Z) dense is uniquely determined by the high density behaviour of the nuclear 


symmetry energy |lf|. Therefore, the 7r /tt + ratio can be used to probe sensitively the EOS 
of neutron-rich matter. On the other hand, within the statistical model for pion production 


2TH , the 7 r /ir + ratio is proportional to exp [(/i n — p p )/T], where T is the temperature, p r 


and [i p are the chemical potentials of neutrons and protons, respectively. At modestly high 
temperatures (T > 4 MeV), the difference in the neutron and proton chemical potentials 
can be written as m 




A P = V asy ~ V asy ~ ^Coulomb + T 


Pn , ^ m + 1 A t 




) m {Pn~PZ) 


m mV 2 ' yrn rp 


( 2 ) 


where Vcouiomb is the Coulomb potential for protons, A t is the thermal wavelength of a 
nucleon and b' s are the inversion coefficients of the Fermi distribution function IE3. The 


difference in neutron and proton symmetry potentials V™ sy — V p sy = 2 v asy (p)S, where the 
function v asy (p) is completely determined by the density-dependence of the symmetry energy 
T5| . It is seen that the kinetic part of the difference p n — p p relates directly to the isospin 


asymmetry p n /p p or p n — p p . Thus within the statistical model too, the /n + ratio is 
sensitive to the ( N/Z)d enae . Moreover, the value of 7r _ /7r + ratio is affected by the competition 
of the symmetry and Coulomb potentials which all depend on the isospin asymmetry of the 
reaction system. These expectations based on two idealized models illustrate qualitatively 
the usefulness of the n~/ t[ + ratio for investigating the EOS of neutron-rich matter. For 
more quantitative studies, however, advanced transport model calculations are necessary. 

To investigate quantitatively effects of the isospin asymmetry of the reaction system on 
the 7T — /7r + ratio, we extend the study of ref. |15| based on a hadronic transport model |22 


We compare three reaction systems 112 Sn + 112 Sn, 124 Sn + 124 Sn and 132 Sn + 124 Sn in the 
beam energy range of 200 to 2000 MeV/nucleon. Since we are interested in situations where 
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the average N/Z ratio of the participant region is the same as that of the initial reaction 
system, we concentrate on central collisions at an impact parameter of 1 fm. To explore the 
isospin-dependence of the nuclear EOS, two phenomenological equations of states with the 
same symmetric part e(p, 0) corresponding to a compressibility Kq = 210 MeV but different 
forms for the symmetry energy E sym (p ) are used. Above the normal nuclear matter density 
Po, the density dependence of E sym (p ) is very controversial. Theoretical results with various 
models and interactions can be classified into two approximately equally large groups, i.e., 
a group where the E sym (p) rises monotonously and one in which it begins to fall above 
about 2p 0 , see e.g., |T3|j2^ , [25|] . As in ref. H, we use the following two representative 
parameterizations 


E S y m(p) — E sym (p 0 )u 


(3) 


and 


sym 


(p) — E sym (p 0 )u 


3 — u 


(4) 


where u = p/p 0 is the reduced density. The linearly increasing E* (p) is the typical 
prediction of the Relativistic Mean Field (RMF) theory |p~4],[26| . The E b sym (p ) first increases 
with density, then starts to decrease above 1.5po and becomes zero at 3po- It approximates 
the prediction of Hartree-Fock calculations using several Skyrme and/or Gogny effective 
interactions [[L3|,|4|,|5|] . 

Shown in Fig. 1 are typical results for the three reactions. In the left panels we examine 
the average isospin asymmetry ( n/p) p > Po of the whole space where the local densities are 
higher than the po- Effects of the symmetry energy is clearly seen after about 10 fm/c when 
a sufficiently high compression has been reached. The E a syrn (p) (E b ym (p)) makes the high 
density region more neutron poor (rich) compared to the initial system due to the continuous 
isospin fractionation |L5|. This effect is stronger for the more neutron-rich systems as one 
expects. Shown in the right panels are the {ji~/'K + )uke ratio 

7 T~ + A' 


(E /tT )/ifce — 


1A° 


7T+ + A ++ + |A+ ’ 


(5) 
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This ratio naturally becomes the final ratio at the freeze-out after all A resonances 

have decayed. In the early stage of the reaction, the (n~/n + )uk e ratio is rather high for the 
neutron richer systems because of the large numbers of neutron-neutron scatterings near 
the surfaces where the neutron skins of the colliding nuclei overlap. It saturates after about 
25 frn/c for all three systems. The higher ir~/ir + ratio with the symmetry energy E b sym (p ) 
reflects directly the higher n/p ratio reached in the high density region. To further examine 
the system dependence, we plot the saturated final n~/i: + ratio as a function (N/Z) system 
in the left window of Fig. 2. As references, we also plotted the ( N/Z) and ( NjZ) 2 . First 
of all, it is seen that the ir~/n + ratio falls far below the first-chance A resonance model 
prediction (N/Z) 2 . This is because of the pion reabsorptions and rescatterings (7r + N A 
and N + A <-► N + N) which reduce the sensitivity of the n~/7T + ratio to the (N/Z) system . 
Moreover, what is more important for the n~/n + ratio is the local, changing n/p ratio during 
particularly the compression phase of the reaction. Secondly, the effect of the symmetry 
energy on the n~/n + ratio is seen to increase slightly as one goes from 112 Sn + 112 Sn to 
124 Sn+ 124 Sn. As one goes further to the 132 SVi+ 124 Sn system, the effect remains at about 
15%. Thus, as far as the n~/ ti + probe of the symmetry energy is concerned, neutron-rich 
stable beams, such as 124 Sn, seem to be sufficient. It is thus tempting to suggest a re¬ 
analysis of existing n~/ tc + data from heavy-ion collisions at GSI energies. The beam energy 
dependence of the n~/n + ratio for the reaction of 132 Sn+ 124 Sn is shown in the right window 
of Fig. 2. While the tt~/ti + ratio decreases with the increasing beam energy, its sensitivity to 
the symmetry energy remains about the same. Similar results are found also for the other 
two reaction systems. The decreasing n~/n + ratio is mainly because of the increasingly 
important contributions of pions from second-chance nucleon-nucleon collisions as the beam 
energy increases. If a first chance nucleon-nucleon collision converts a neutron to a proton 
by producing a n ~, subsequent collisions of the still energetic proton can convert itself back 
to a neutron by producing a n + . Eventually, at very high energies the sequential multiple 
nucleon-nucleon collisions will lead to a n~/rc + « 1. 

In summary, within an isospin-dependent hadronic transport model using two repre- 
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sentative density-dependent symmetry energy functions predicted by many-body theories, 
it is shown that the n~/n + ratio increases with the isospin asymmetry, but decrease with 
the beam energy of the reaction. For neutron-rich systems, such as, 124 Sn + 124 Sn and 
132 Sn + 124 Sn, the sensitivity of n~/n + ratio to nuclear symmetry energy is almost beam 
energy independent. 

The author would like to thank W.G. Lynch for suggesting me to conduct this research 
and useful discussions. This work was supported in part by the National Science Foundation 
Grant No. PHY-0088934 and Arkansas Science and Technology Authority Grant No. 00-B- 
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FIGURES 



FIG. 1. Left panels: average neutron/proton ratio in the whole space with densities higher 
than the normal nuclear matter density as a function of time. Right panels: the {pi~/TT + )uke ratio 
as a function of time for the three reactions. The solid (dashed) lines are the results using the 
symmetry energy E“ ym ( P ) (. E b sym {p )) 
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FIG. 2. The (pi /vr + ) ratio as a function of the isospin asymmetry (left window) and beam 
energy (right window) of the reaction system 
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